Noise of a Quantum-Dot System in the Cotunneling Regime 



O 

o 

C 

a 

Hi 

(N 
(N 



Eugene V. Sukhorukov, Guido Burkard, and Daniel Loss 
Department of Physics and Astronomy, University of Basel, 
Klingelbergstrasse 82, CH-4056 Basel, Switzerland 

We study the noise of the cotunneling current through one or several tunnel-coupled quantum dots 
in the Coulomb blockade regime. The various regimes of weak and strong, elastic and inelastic 
cotunneling are analyzed for quantum-dot systems (QDS) with few-level, nearly-degenerate, and 
continuous electronic spectra. We find that in contrast to sequential tunneling where the noise is 
either Poissonian (due to uncorrelated tunneling events) or sub-Poissonian (suppressed by charge 
conservation on the QDS) , the noise in inelastic cotunneling can be super-Poissonian due to switching 
between QDS states carrying currents of different strengths. In the case of weak cotunneling we 
prove a non-equilibrium fluctuation-dissipation theorem which leads to a universal expression for 
the noise-to-current ratio (Fano factor). In order to investigate strong cotunneling we develop a 
microscopic theory of cotunneling based on the density-operator formalism and using the projection 
operator technique. The master equation for the QDS and the expressions for current and noise 
in cotunneling in terms of the stationary state of the QDS are derived and applied to QDS with a 
nearly degenerate and continuous spectrum. 

PACS numbers: 73.23.-b, 73.23.Hk, 72.70.+m, 73.63.Kv, 73.63.-b 
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I. INTRODUCTION 



In recent years, there has been great interest in 
transport properties of strongly interacting mesoscopic 
systemsJd As a rule, the electron interaction effects be- 
come stronger with the reduction of the system size, since 
the interacting electrons have a smaller chance to avoid 
each other. Thus it is not surprising that an ultrasmall 
quantum dot connected to leads in the transport regime, 
being under additional control by metallic gates, pro- 
vides a unique possibility to study strong coraelation ef- 
fects both in the leads and in the dot itself.0 This has 
led to a large number of publications on quantum dots, 
which investigate situations where the current acts as a 
probe of correlation effects. Historically, the nonequilib- 
rium current fluctuations (shot noise) were initially con- 
sidered as a serious problem for device applications of 
quantum dotsu El rather than as a fundamental physical 
phenomenon. Later it became clear that shot noise is an 
interesting phenomenon in itself p because it contains ad- 
ditional information about correlations, which is not con- 
tained, e.g., in the linear response conductance and can 
be used as a further approach-tQ study transport in quatu 
turn dots, both theoreticallyQd'Ql£3 and experimentally.^ 

Similarly, the majority of papers on the noise of quan- 
tum dots consider the sequential (single-electron) tunnel- 
ing regime, where a classical description (the so-called 
"orthodox" theory) is applicable. cfl We are not aware 
of any discussion in the literature of the shot noise in- 
duced by-ira-i cotunneling (two-electron, or second-order) 
current 0tB except Ref. ^l], where the particular case 
of weak cotunneling (see below) through a double-dot 
(DD) system is considered. Again, this might be because 
until very recently cotunneling has been regarded as a 
minor contribution to the sequential tunneling current, 



which spoils, the precision of single-electron devices due 
to leakageEH However, it is now well understood that co- 
tunneling is interesting in itself, since it is responsible 
for strongly corr elate d effects such as the Kondo effect 
in quantum dotsj23'E2l or can be used as, a probe of two- 
electron entanglement and nonlocalityj£il etc. 

In this paper we present a thorough analysis of the 
shot noise in the cotunneling regime. Since the single- 
electron "orthodox" theory cannot be applied to this 
case, we first develop a microscopic theory of cotunneling 



suitable for the calculation of the shot noise in Sees. Ill 



and IV. [For an earlier microscopic theory of transport 
through quantum dots see Refs. |30|-p^.] We consider the 



transport through a quantum-dot system (QDS) in the 
Coulomb blockade (CB) regime, in which the quantiza- 
tion of charge on the QDS leads to a suppression of the 
sequential tunneling current except under certain reso- 
nant conditions. We consider the transport away from 
these resonances and study the next-order contribution 
to the current, the so-called cotunneling current P 5 ^ 6 ! In 
general, the QDS can contain several dots, which can be 
coupled by tunnel junctions, the double dot (DD) being a 
particular examplcO The QDS is assumed to be weakly 
coupled to external metallic leads which are kept at equi- 
librium with their associated reservoirs at the chemical 
potentials /i;, I = 1,2, where the currents // can be mea- 
sured and the average current / through the QDS is de- 
fined by Eq. (EDI). 



Before proceeding with our analysis we briefly review 
the results available in the literature on noise of sequen- 
tial tunneling. For doing this, we introduce right from 
the beginning all relevant physical parameters, namely 
the bath temperature T, bias A/i = fix — [1%, charg- 
ing energy Ec, average level spacing 8E, and the level 
width r — Ti +T2 of the QDS, where the tunneling rates 
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FIG. 1. Schematic representation of the quantum dot sys- 
tem (QDS) coupled to two external leads 1 and 2 (light grey) 
via tunneling barriers (dark grey), where the energy scale is 
drawn vertically. The tunneling between the QDS and the 
leads I = 1,2 is parametrized by the tunneling amplitudes 
71, where the lead and QDS quantum numbers k and p have 
been dropped for simplicity, see Eq. (2.2). The leads are at the 
chemical potentials /Ui,2, with an applied bias A/j, = /xi — /xz- 
The (many-particle) eigenstates of the QDS with one added 
electron (N + 1 electrons in total) are indicated by their ener- 
gies Ei, E2, etc., with average level-spacing SE. The energy 
cost for adding a particle from the Fermi level of lead I to the 
N-electron QDS is denoted by A+(Z, N) > and is strictly 
positive in the CB regime. Note that the energies A-(l,N) 
for removing particles from the QDS containing N electrons 
are positive as well, and are not drawn here. The cotunneling 
process is visualized by two arrows, leading from the initial 
state in, say, lead 1 (full circle), via a virtual state on the 
QDS (open circle), to the final state in lead 2 (full circle). 



Ti = ttv\Ti\ 2 to the leads 1 — 1,2 are expressed in terms of 
tunneling amplitudes T) and the density of states u eval- 
uated at the Fermi energy of the leads. In Fig. |l| the most 
important parameters are shown schematically. This va- 
riety of parameters shows that many different regimes 
of the CB are possible. In the. linear response regime, 
A/i <C ksT, the thermal noiseH is given by .-the equilib- 
rium fluctuation-dissipation theorem (FDT)Ea Although 
the cross-over from the thermal to noncquilibrium noise 
is of our interest (see Sec. [II ), in this section we discuss 
the shot noise alone and set T = 0. Then the noise at 
zero frequency to = 0, when 5I 2 = —SI\, can be charac- 
terized by one single parameter, the dimensionless Fano 
factor F = S(0)/e\I\, where the spe ctra l density of the 
noise S(0) = S 22 (0) is defined by Eq. The Fano fac- 

tor acquires the value F = 1 for uncorrclatcd Poissonian 
noise. 

Next we discuss the different CB regimes. (1) In the 
limit of large bias A/x ^> Eq, when the CB is sup- 
pressed, the QDS can be viewed as being composed of 
two tunnel junctions in series, with the total conduc- 



tance G = GiG 2 /(Gi + G 2 ), where G = 7re 2 w D |T ; | 2 is 
the conductance of the tunnel junctions to lead I, and 
vd is the density of dot states. Then the Fano fac- 
tor is given by F = (G? + G\)/(Gi + G 2 ) 2 , as it has 
been found in Refs. f|,||,[7| Thus, the shot noise is sup- 
pressed, F < 1, and reaches its minimum value for the 
symmetric QDS, G x = G 2 , where F = 1/2. (2) The 
low bias regime, SE <C A/x <C Ec- The first inequal- 
ity SE <C A/x allows to assume a continuous spectrum 
on of the QDS and guarantees that the single-electron 
"orthodox" theory based on a classical master equation 
can be applied. The second inequality A/x <C Ec means 
that the QDS is in the CB regime, where the energy cost 
A±(i, iV) = E(N ± 1) - E(N) T Mi f or the electron tun- 
neling from the Fermi level of the lead I to the QDS (+) 
and vice versa (— ) oscillates as a function of gate voltage 
between its minimum value A± < (where the energy 
deficit turns into a gain, |A±| ~ A/x) and its maximum 
value A± ~ Ec- Here, E(N) denotes the ground-state 
energy of the A-electron QDS. Thus the current I as 
a function of the gate voltage consists of the CB peaks 
which are at the degeneracy points A± < 0, where the 
number of electrons on the QDS fluctuates between N 
and N + 1 due to single-electron tunneling. The peaks 
are separated by plateaus, where the single-electron tun- 
neling is blocked because of the finite energy cost A± > 
and thus the sequential tunneling current vanishes. At 
the peaks the current is given by I — rf&iffe / (71 +72), 
while the Fano factor has been reportecHlJ~LL3 to be equal 
to F = (7? + 7 2 2 )/(7i + 7 2 ) 2 , 1/2 < F < 1, where 
71 = e- 2 Gi|A+(l,A)| and l2 = e~ 2 G 2 |A_(2, A + 1)| 
are the tunneling rates to the QDS from lead 1 and from 
the QDS to lead 2, respectively. Within the "orthodox" 
theory tunneling is still possible between the peaks at 
finite temperature due to thermal activation processes, 
and then the Fano factor approaches the Poissonian value 
F = 1 from below. (3) Finally, the limit r <C A/i SE 
is similar to the previous case, with the only difference 
that the dot spectrum is discrete. The sequential tunnel- 
ing picture can still be applied; the result for the Fano 
factor at the current peak is. F = (Tj + r^)/(r x + T 2 ) 2 , 
so that again 1/2 < F < 1.E3 

We would like to emphasize the striking similarity of 
the Fano factors in all three regimes, where they also 
resemble the Fano factor of the noninteracting double- 
barrier system. U The Fano factors in the first and second 
regimes become even equal if the ground-state level of 
the QDS lies exactly in the middle between the Fermi 
levels of lead 1 and 2, |A+| = |A_|. We believe that 
this "ubiquitous"Q double-barrier character of the Fano 
factor can be interpreted as being the result of the natu- 
ral correlations imposed by charge conservation rather 
than by interaction effects. Indeed, in the transport 
through a double-barrier tunnel junction each barrier can 
be thought of as an independent source of Poissonian 
noise. And although in the second regime the CB is ex- 
plicitly taken into account, the stronger requirement of 
charge conservation at zero frequency, SI\ + SI 2 = 0, has 
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to be satisfied, which leads to additional correlations be- 
tween the two sources of noise and to a suppression of the 
noise below the Poissonian value. At finite frequency (but 
still in the classical range defined as w C A/j,,Ec) tem- 
porary charge accumulation on the QDS is allowed, and 
for frequencies larger than the tunneling rate, lo 3> 71,2, 
the conservation of charge does not need to be satisfied, 
while the noise power S22 approaches its Poissonian valua 
from below, and the cross correlations vanish, S12 = OSB 
Based on this observation we expect that the direct mea- 
surement of interaction effects in noise is onljf possible 
either in the quantum (coherent) CB regimeta A/i ~ T 
or in the Kondo regime,EZni3 where both charge conser- 
vation and many-electron effects lead to a suppression of 
the noise. Another example is the noise in the quantum 
regime, A/i < lo ~ Ec, where it contains singularities as- 
sociated with the . "aho ta- assisted transitions" above the 
Coulomb gap A±.eMeI 

To conclude our brief review we would like to em- 
phasize again that while the zero-frequency shot noise 
in the sequential tunneling regime is always suppressed 
below its full Poissonian value as a result of charge con- 
servation (interactions suppressing it further), we find 
that, in the present work the shot noise in the cotun- 
neling regimetll is either Poissonian F — 1 (elastic or 
weak inelastic cotunneling) or, rather surprisingly, non- 
Poissonian F ^ 1 (strong inelastic cotunneling). There- 
fore the non-Poissonian noise in QDS can be considered 
as being a fingerprint of inelastic cotunneling. This dif- 
ference of course stems from the different physical origin 
of the noise in the cotunneling regime, which we dis- 
cuss next. Away from the sequential tunneling peaks, 
A± > 0, single-electron tunneling is blocked, and the 
only elementary tunneling process which is compatible 
with energy conservation is the siipwlteneous tunneling 
of two electrons called cotunnelingE3c3. In this process 
one electron tunnels, say, from lead 1 into the QDS, and 
the other electron tunnels from the QDS into lead 2 with 
a time delay on the order of Ajjj 1 (see Ref. ^l|). This 
means that in the range of frequencies, lo <C A±, (which 
we assume in our paper) the charge on the QDS does not 
fluctuate, and thus in contrast to the sequential tunneling 
the correlation imposed by charge conservation is not rel- 
evant for cotunneling. Furthermore, in the case of elastic 
cotunneling (A/i < 6E), where the state of the QDS re- 
mains unchanged, the QDS can be effectively regarded as 
a single barrier. Therefore, subsequent elastic cotunnel- 
ing events are uncorrelated, and the noise is Poissonian 
with F = 1. On the other hand, this is not so for inelas- 
tic cotunneling (A/i > 6E), where the internal state of 
the QDS is changed, thereby changing the conditions for 
the subsequent cotunneling event. Thus, in this case the 
QDS switches between different current states, and this 
creates a correction to noise AS, so that the total noise is 
non-Poissonian, and can become super-Poissonian. The 
other mechanism underlying super-Poissonian noise is 
the excitation of high energy levels (heating) of the QDS 
caused by multiple inelastic cotunneling transitions and 



leading to the additional noise ASh- Thus the total noise 
can be written as S = el + ASh + AS. For other cases 
exhibiting super-Poissonian noise (in the strongly non- 
linear bias regime) see Ref. @. 

According to this picture we consider the following 
different regimes of the inelastic cotunneling. We first 
discuss the weak cotunneling regime w <C Win, where 
W ~ riF2A/i/Aj_ is the average rate of the inelastic co- 



tunneling transitions on the QDS [see Eqs. (4.23- 4~26| )], 
and Wi n is the intrinsic relaxation rate of the QDS to its 
equilibrium state due to the coupling to the environment. 
In this regime the cotunneling happens so rarely that the 
QDS always relaxes to its equilibrium state before the 
next electron passes through it. Thus we expect no cor- 
relations between cotunneling events in this regime, and 
the zero-frequency noise is going to take on its Poisso- 
nian value with Fano factor F — 1, as first obtained for 
a special case in Ref. ^l|. This result is generalized in 
Sec. Ill, where wc find a universal relation between noise 



and current of single-barrier tunnel junctions and, more 
generally, of the QDS in the first nonvanishing order in 
the tunneling perturbation V. Beca use of the universal 
character of the results Eqs. (3.10) and ( |3.21 ) we call 
them the nonequilibrium FDT in analogy with linear re- 
sponse theory. 

Next, we consider strong cotunneling, i.e. w 3> wi n . 
The microscopic theory of the transport and noise in 
this regime based on a projector operator technique is 
develop ecLiiL Sec. IV. In the case of a few-level QDS, 
SE ~ Ecf^i noise turns out to be non-Poissonian, as we 
have discussed above, and this effect can be estimated as 
follows. The QDS is switching between states with the 
different currents I ~ ew, and we find SI ~ ew. The 



QDS stays in each state for the time r 



There- 



fore, for the positive correction to the noise power we get 
AS ~ SI 2 t ~ e 2 w, and the estimate for the correction to 
the Fano factor follows as AS /el ~ 1. A similar result 
is expected for the noise induced by heating, ASh, which 
can roughly be estimated by assuming an equilibrium dis- 
tribution on the QDS with the temperature k B T ~ Au. 
and considering the additional noise as being thermalJEI 
AS h ~ Gk B T ~ (eI/A(i)k B T ~ el. The characteris- 
tic frequency of the noise correction AS is lo ~ w, with 
AS vanishing for lo 3> w (but still in the classical range, 
lo <C A/i). In contrast to this, the additional noise due 
to heating, ASh, does not depend on the frequency. 

In Sec. |v| we consider the particular case of nearly de- 
generate dot states, in which only few levels with an en- 
ergy distance smaller than SE participate in transport, 
and thus heating on the QDS can be neglected. Specif- 
ically, for a two-level QDS we predict giant (divergent) 
super-Poissonian noise if the off-diagonal transition rates 
vanish. The QDS goes into an unstable mode where it 
switches between states 1 and 2 with (generally) different 
currents. We consider the transport through a double- 
dot (DD) system as an example to illustrate this effect 
[see Eq. ( pT2| ) and Fig. §. 

Finally, we discuss the case of a multi-level QDS, 
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SE <C Eq- In this case the correlations in the cotun- 
neling current described above do not play an essential 
role. In the regime of low bias, Ap <^MLE Ec) 1 ^ 2 , elas- 
tic cotunneling dominates transport ,E£Il3 and thus the 
noise is Poissonian. In the opposite case of large bias, 
A± ^> A/i ^> (SE Ec) 1 / 2 , the transport is governed by 
inelastic cotunneling, and in Sec. VI we study heating 
effects which are relevant in this regime. For this we use 
the results of Sec. IV and derive a kinetic equation for 



the distribution function /(e). We find three universal 
regimes where I ~ A/j, 3 , and the Fano factor does not 
depend on bias the A/i. The first is the regime of weak 
cotunneling, r- m <C r c , where r- m and r c are time scales 
characterizing the single-particle dynamics of the QDS. 
The energy relaxation time T m describes the strength of 
the coupling to the environment while r c ~ eisoAp/I 
is the cotunneling transition time. Then we obtain for 
the distribution f(s) = 9(—e), reproducing the result of 
Ref. We also find that F = 1, in agreement with the 
FDT proven in Sec. |l|. The other two regimes of strong 
cotunneling r m t c are determined by the electron- 
electron scattering time r ee . For the cold-electron regime, 
t c <C T ee , we find the distribution function by solving the 
integral equations (6.11) and ( |6.12 ), while for hot elec- 
trons, t c 3> T ee , / is given by the Fermi distribution 
function with an electron temp erature obtained from the 
energy balance equation ( |6.15| ). We use f(e) to calculate 
the Fano factor, which turns out to be very close to 1. On 
the other hand, the current depends not only on G\Gi 
but also on the ratio, G1/G2, depending on the cotun- 
neling regime [see Fig. 4]. Details of the calculations are 
deferred to four appendices. 



II. MODEL SYSTEM 

The quantum-dot system (QDS) under study is weakly 
coupled to two external metallic leads which are kept in 
equilibrium with their associated reservoirs at the chem- 
ical potentials //;, I — 1,2, where the currents can 
be measured. Using a standard tunneling Hamiltonian 
approach,c2l we write 



H = H Q + V, 
Hr = W 



1 = 1,2 k 



Hq = Hl + Hs + H[ nt 

t - 



/ Spdpdp . 
p 



V = (A + d\) , A = E r '*p 4 d p . 



(2.1) 
(2.2) 

(2.3) 



1=1.2 



k , p 



where the terms Hl and Hs describe the leads and QDS, 
respectively (with k and p from a complete set of quan- 
tum numbers), and tunneling between leads and QDS is 
described by the perturbation V. The interaction term 
Hi n t is specified below. The iV-electron QDS is in the 
cotunneling regime where there is a finite energy cost 
A±(l,N) > for the electron tunneling from the Fermi 
level of the lead I to the QDS (+) and vice versa (— ), so 
that only processes of second order in V are allowed. 



To describe theptransport through the QDS we apply 
standard methodscJ and adiabatically switch on the per- 
turbation V in the distant past, t = to — > —00. The 
perturbed state of the system is described by the time- 
dependent density matrix p(t) 
which can be written as 



e -iff(f-t )p giff(i-to) ) 



p(t) = e-^-^po , LA = [H, A], MA , (2.4) 

with the help of the Liouville operator L = Lq + 
Here po is the grand canonical density matrix of the un- 
perturbed system, 



Po 



z -X e -K/k B T 



(2.5) 



where we set K = H — J2i V1N1. 

Because of tunneling the total number of electrons in 
each lead Ni = ^2 k cLcjfe is no longer conserved. For the 
outgoing currents // = eNi we have 

It =ei[V,Ni\ =ei(D\- D t ) . (2.6) 

The observables of interest are the average current / = 
I2 = —Ii through the QDS, and the spectral density of 
the noise Sw(u!) = J dtSw(t) exp(iu>t), 

h = Trp(0)i ; , S w (t) = ReTr p(0)SIi(t)SIv(0) , (2.7) 

where SIi = J; — Jj. Bel ow we will use the interaction rep- 
resentation where Eq. ( |2.7|) can be rewritten by replacing 
p(0) -> po and /,(t) W(t)Ii(t)U(t), with 



U(t) = Texp 



dt' V(t') 



(2. 



In this representation, the time dependence of all opera- 
tors is governed by the unperturbed Hamiltonian Ho- 



III. NON-EQUILIBRIUM 
FLUCTUATION-DISSIPATION THEOREM FOR 
TUNNEL JUNCTIONS 

In this section we prove the universality of noise of tun- 
nel junctions in the weak cotunneling regime w <C uii n 
keeping the first nonvanishing order in the tunneling 
Hamiltonian V. Since our final results Eqs. ( ft.lOp , ( |3.12fl , 
(fyH ), and ( |3~2l| ) can be applied to quite general systems 
out-of-equilibrium we call this result the non-equilibrium 
fluctuation-dissipation theorem (FDT). In particular, the 
geometry of the QDS and the interaction Hi n t are com- 
pletely arbitrary for the discussion of the non-equilibrium 
FDT in this section. Such a non-equilibrium-J^DT was 
derived for single barrier junctions long ago.c3 We will 
need to briefly review this case which allows us then to 
generalize the FDT to QDS considered here in the most 
direct way. 
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A. Single-barrier junction 

The total Hamiltonian of the jun ctio n [given by 

Eqs. and the currents Eq. have to be 

replaced by H = Hl + H lnt + V, where 



V = A + A*, A = J2 T kk'4 k ci k> 



k,k' 



h = -h = ei [V 7 N 2 ] = ei (A* - A) 



(3.1) 



(3.2) 



For the sake of generality, we do not specify the interac- 
tion _ff; nt in this section, nor the electron spectrum in the 
leads, and the geometry of our system. 

Applying-, the standard interaction representation 



technique)!!] we expand the expression (2.8) for U(t) and 
keep only first non-vanishing contributions in V , obtain- 



ing 



I(t) 



dt'i 



V(t'),h(t) ) 



(3.3) 



where we use the notation (...) = Trpo(. ..). Analo- 
gously, we find that the first non-vanishing contribution 
to the noise power S(u>) = S 22 (lo) is given by 



dte iut ({I 2 (t)J 2 (0)}) 



(3.4) 



I = e I dt([Ai(t),A(p)]), 

— OO 

OO 

S{u) = e 2 { dtcos(ujt)({A^(t),A(0)}) 



where {. . .} stands for anticommutator, and Jf = in 
leading order. 

We notice that in Eqs. (O) and (O) the terms (AA) 
and (A^ A^) are responsible for Cooper pair tunneling and 
vanish in the case of normal (interacting) lead s. Taking 
this into account and using Eqs. (pj) and (|3.2[) we obtain 



(3.5) 
(3.6) 



where we also used (At(t)A(Q)) = (A^ (O)A(-t)). 

Next we app ly t he spec tral decomposition to the cor- 
relators Eqs. (3.5) and (|3.6| ), a similar procedure to 
that which also leads to the equilibrium fluctuation- 
dissipation theorem. The crucial observation is that 
[H ,Ni] = 0, I = 1,2 (we stress that it is only 
the tunneling Hamiltonian V which does not commute 
with Ni, while all interactions do not change the num- 
ber of electrons in the leads). Therefore, we are al- 
lowed to use for our spectral decomposition the ba- 
sis |n) = \E n ,Ni,N 2 ) of eigenstates of the operator 
K = Hq — X^jU/iV/, which also diagonalizes the grand- 
canonical density matrix po [given by Eq. (|2.5|)], p n — 



(n|po| n ) = Z 1 cxp[— E n /k B T]. Next we introduce the 
spectral function, 



A{oj) = 27r^(p n + p m )|(m|A|n) 



x5(lv + E n 



Er, 



(3.7) 



and rewrite Eqs. (3.5) and ( |3.6[ ) in the matrix form in the 
basis |n) taking into account that the operator A c reat es 
(annihilates) an electron in the lead 2 (1) [see Eq. (3.1)]. 
We obtain following expressions 



e tanh 

,2 



A/i 



2k B T 



A(Afi) . 



(3.8) 
(3.9) 



where A/i = fi\ — /i 2 . From these equations our main 
result follows 



S(u,Afi) = |^coth 



A/i ± U) 



2k B T 



I(Ah±uj), (3.10) 



where we have neglected contributions of order 
An/EF,Qj/eF <C 1. We call the relation ( |3 .10 ) non- 
equilibrium fluctuation-dissipation theorem because of 
its general validity (we recall that no assumptions on ge- 
ometry or interactions were made). 

The fact that the spectral function Eq. ( |3.7| ) depends 
only on one parameter can be used to obtain further 
useful relations. Suppose that in addition to the bias 
A/i a small perturbation of the form 5/ie~' luJt is applied 
to the junction. This perturbation generates an ac cur- 
rent 6I(lo, A/i)e _iwt through the barrier, which depends 
on both parameters, to and A/i. The quantity of in- 
terest is the linear response conductance G(to, A/i) = 
eSI(Lu, A/i)/<5/i. The perturbation Sfi can be taken into 
account in a standard way by multiplying the tunnel- 
ing amplitude A(t) by a phase factor e~ l ^\ where 6 = 
S/ie~ luJt . Substituting the new amplitude into Eq. (3.3) 



and expanding the current with respect to Sfi, we arrive 
at the following result, 



ReG(w, A/t) 



dtsm{ut){[A\t),A{Q)]) . (3.11) 



Finally, applying the spectral decomposition to this equa- 
tion we obtain 



(2/e)w Re G(w, A/i) = /(A/i 



/(A/i 



(3.12) 



which holds for a general nonline ar / v s A/i dependence. 
From this equation and from Eq. (3.10) it follows that the 
noise power at zero frequency can be expressed through 
the conductance at finite frequency as follows 



5(0, A/i) + 5(0, -A/i) 



2lo coth 



2k B T 



ReG(w,0)|^ AM . (3.13) 
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And for the noise power at zero bias we obtain S(u>, 0) = 
u> coth(gj/2fcBy4Re G(u , 0), which is the standard equi- 
librium FDTHj Eq. ( |3 .10 ) reproduces the result of 
Ref. The current is not necessary linear in A/j, (the 
case of tunneling into a Luttingcr liquidc3 is an obvi- 
ous example) , and in the limit T, uj — > we find the 
Poissonian noise, 5* = el. In the limit T, A/i — > 0, the 
quantum noise becomes S(oj) — e[I(ui) — I(—ui)]/2. If 
/( — A/i) = —I(A[i), we get S(ui) = e/(w), and thus S(cu) 
can be obtained from /(A/i — > uj). 



B. Quantum dot system 

We consider now tunneling through a QDS. In this 
case the problem is more complicated: In general, the 
two currents Ii are not independent, because [h^h] ^ 0, 
and thus all correlators Sw are nontrivial. In particular, 
it has been proven in Ref. 21 that the cross-correlations 
ImS'i2(w) are sharply peaked at the frequencies uo — A±, 
which is caused by a virtual charge-imbalance on the 
QDS during the cotunncling process. The charge accu- 
mulation on the QDS for a time of order Ajjj 1 leads to an 



additional contribution to the noise at finite frequency uj. 
Thus, we expect that for uj ~ A± the correlators Sw can- 
not be expressed through the steady-state current I only 
and thus / has to be complemented by some other dis- 
sipative counterparts, such as differential conductances 
G w (see Sec. pIA| ). 

On the other hand, at low enough frequency, to <C A±, 
the charge conservation on the QDS requires SI S — 
(SI2 + 5I\)/2 sa 0. Below we concentrate on the limit 
of low frequency and neglect contributions of order of 
oj/A± to the noise power. In Appendix [A| we prove that 
S ss ~ (lu/A±) 2 , and this allows us to redefine the cur- 
rent and the noise power as I = Id = {I2 — -M/2 and 
S(lu) = 5dd(w).0 In addition we require that the QDS 
is in the cotunneling regime, i.e. the temperature is low 
enough, fc^T <C A±, although the bias A/j, is arbitrary 
(i.e. it can be of the order of the energy cost) as soon as 
the sequential tunneling to the dot is forbidden, A± > 0. 
In this limit the current through a QDS arises due to the 
direct hopping of an electron from one lead to another 
(through a virtual state on the dot) with an amplitude 
which depends on the energy cost A± of a virtual state. 
Although this process can change the state of the QDS, 
the fast energy relaxation in the weak cotunneling regime, 
w <C uhn, immediately returns it to the equilibrium state 
(for the opposite case, see Sees. IV- VI). This allows us to 
apply a perturbation expansion with respect to tunnel- 
ing V and to keep only first nonvanishing contributions, 
which we do next. 

It is convenient to introduce the notation Di(t) = 
J_ dt' Di(t'). We notice that all relevant matrix ele- 
ments, (N\Di(t)\N + 1) ~ e- lA +*, (N - l\Di(t)\N) - 
e , are fast oscillating functions of time. Thus, under 
the above conditions we can write D;(oo) = 0, and even 



more general, J^°° dt Di(t)e ±tuJt = (note that we have 
assumed earlier that uj -C A±). Using these equalities 
and the cyclic property of the trace we obtain the follow- 
ing result (for details of the derivation, see Appendix [A|) , 



I = e / dt{[B\t),B{0)]), 



B = D 2 D\ + d\d 2 . 



(3.14) 
(3.15) 



Applying a similar procedure (see Appendix A), we arrive 
at the following expression for the noise power S = S22, 
see Eq. (O), 



S(w) 



dt cos(ujt)({B\t),B(0)}) . (3.16) 



where we have dropped a small contribution of order 
w/A±. 



Thus, we have arrived at Eqs. ( 3.14 ) a nd (3.16 ) which 
are formally equivalent to Eqs. (3.5) and (3.6). Similarly 
to A in the single-barrier case, the operator B plays the 
role of the effective tunneling amplitude, which annihi- 
lates an electron in lead 1 and creates it in lead 2. Similar 



to Eqs. (3/7), (3^8), and (3J5) we can express the current 



and the noise power 



I(Afi) — e tanh 



Afi 



2k B T 



£(A M ) , 



S( W ,A M ) = y^B(A/i±a;) ! 



(3.17) 
(3.18) 



in terms of the spectral function 

6H = 27r^(p n + p m )|<m|B|n)| 2 

n,m 

xSiiU + En-Ern). (3.19) 

The difference, however, becomes obvious if we notice 



that in contrast to the operator A [see Eq. (3.1)] which is 
a product of two fermionic Schrodinger operators with an 
equilibrium spectrum, the operator B contains an addi- 
tional time integration with the time evolution governed 
by Ho = K + J^i HiNi- Applying a further sp ectra l de- 
composition to the operator B [given by Eq. ( 3.15| )1 we 
arrive at the expression 



'H%)=ET 2 ' n l (n ^ |n 

n' 

+ E 



En' — E n — fli 

m|D||n")(n"|£>2|n) 



/'2 



(3.20) 



where the two sums over n' and n" on the Ihs are different 
by the order of tunneling sequence in the cotunneling pro- 
cess. Thus we see that the current and the noise power 
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depend on both chemical potentials pi.2 separately (in 
contrast to the one-para mete r dependence for a single- 
barrier junction, see Sec. Ill A| ) , and therefore the shift of 
A/it in Eq. ( 3.18 ) by ±ui will also shift the energy denom- 
inators of the matrix elements on the Ihs of Eq. (3.20). 
However, since the energy denominators are of order A± 
the last effect can be neglected and we arrive at the final 
result 



S(ui, A/i) = - coth 



Ap ± ui 
2k B T 



I(Ap,±u>) 

+0(w/A±). (3.21) 



This equation represents our nonequilibrium FDT for the 
transport through a QDS in the weak cotunneling regime. 
A special case with T, ui = 0, giving S — el, has been 
derived in Ref. ^l[ To conclude this section we would 
like to list again the conditions used in the deriv ation . 
The universality of noise to current relation Eq. (3.21) 
proven here is valid in the regime in which it is sufficient 
to keep the first nonvanishing order in the tunneling V 
which contributes to transport and noise. This means 
that the QDS is in the weak cotunneling regime with 
ui, ksT <C A±, and w m ^> w. 



IV. MICROSCOPIC THEORY OF STRONG 
COTUNNELING 



A. Formalism 



In this section, we give a systemat ic m icroscopic 
derivation of the m aster equation, Eq. ( 4. 22] ), the av- 
erage curr en t, Eq . (4.37), and the current correlators, 



Eqs. ( |4.52| )-( f4.54| ) for the QDS coupled to leads, as in- 
troduced in Eqs. (2.1)-(2.3), in the strong cotunneling 
regime, Wi n <C w. Under this assumption the intrinsic 
relaxation in the QDS is very slow and will in fact be 
neglected. Thermal equilibrati on ca n only take place via 
coupling to the leads, see Sec. [VB . Due to this slow re- 
laxation in the QDS we find that there arc non-Poissonian 
correlations AS in the current through the QDS because 
the QDS has a "memory" ; the state of the QDS after the 
transmission of one electron influences the transmission 
of the next electron. A basic assumption for the following 
procedure is that the system and bath are coupled only 
weakly and only via the perturbation V, Eq. (pT^). The 
inter acti on part iJi nt of the unperturbed Hamiltonian Hq, 
Eq. (2.1), must therefore be separable into a QDS and a 



lead part, i?; n t = H™* + Hf^. Moreover, Ho conserves 
the number of electrons in the leads, [Ho, Ni] = 0, where 

We assume that in the distant past, to 
system is in an equilibrium state 



-co, the 



Po = PS <X> PL, 



PL 



1 



-K L /k B T 



(4.1) 



where Z L = Tr exp[-K L /k B T], K L = H L - J2iVi N h 
and ui is the chemical potential of lead Note that both 
leads are kept at the same t emp erature T . Physically, the 
product form of po in Eq. (4.1) describes the absence of 
correlations between the QDS and the leads in the initial 
state at to. Furthermore, we assume that the initial state 
Po is diagonal in the eigenbasis of Ho, i.e. that the initial 
state is an incoherent mixture of eigenstates of the free 
Hamiltonian. 

In systems which can be divided into a (small) system 
(like the QDS) and a (possibly large) external "bath" at 
thermal equilibrium (here, the leads coupled to the QDS) 
it turns out to be i vory , u seful to make use of the superop- 
erator formalism, e1Ie3cj and of projectors Pt = PlT^l, 
which project on the "relevant" part of the density ma- 
trix. We obtain PtP by taking the partial trace Tr^ of p 
with respect to the leads and taking the tensor product 
of the resulting reduced density matrix with the equi- 
librium state pl- Here, we will consider the projection 
operators 



P = (PdPn ® Il)Pt, Q = l-P, 



(4.2) 



satisfying P 2 = P, Q 2 = Q, PQ = QP =j-Q, where P is 
composed of Pt and two other projectorsc3 Pjj and Pn, 
where Pd projects on operators diagonal in the eigenba- 
sis {\n)} of Hs, i.e. (n\P£>A\m) — S nm (n\A\m) , and P/v 
projects on the subspace with N particles in the QDS. 
The particle number N is defined by having minimal en- 
ergy in equilibrium (with no applied bias); all other par- 
ticle niupbers have energies larger by at least the energy 
deficito A. Above assumptions about the initial state 
Eq. ( [4.1D of the system at to — > — 00 can now be rewrit- 
ten as 



Ppo = Po- 



(4.3) 



For the purpose of deriving the master equation we 
take the Lap lace transform of the time-dependent density 
matrix Eq. (2.4), with the result 



p(z) = R(z)p . 



(4.4) 



Here, R(z) is the resolvent of the Liouville operator L, 
i.e. the Laplace transform of the propagator exp(— itL), 

R{z)= / dte lt(z ~ L) =i{z-L)- 1 = —!— , (4.5) 
Jo z-L 

where z = u> + ir\. We choose 77 > in order to ensure 
convergence (L has real eigenvalues) and at the end of 
the calculation take the limit 77 — > 0. We can split the 
resolvent into four parts by multiplying it with the unity 
operator P + Q from the left and the right, 



R = PRP + QRQ + PRQ + QRP. 



(4.6) 



Inserting the identity operator —i(z — L)R(z) = —i(z — 
L)(P + Q)R(z) between the two factors on the left hand 
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side of QP = 0, PQ = 0, Q 2 = Q, and P 2 
obtain 



±—QL v PR(z)P, 
z - QLQ 

PR(z)Q = -iPR (z)PL v QR(z)Q, 

i 



QR(z)P = Q 
PR(z)Q = - 
QR(z)Q = Q 



z - QLQ + iQL v PR Q (z)PL v Q 



PR(z)P = P- 



-P, 



z-E(z) 

where we define the self-energy superoperator 

f 



E(z) = PL V Q 



z-QLQ 



QL V P, 



P, we 

(4.7) 
(4.8) 
(4.9) 

(4.10) 
(4.11) 



and the free resolvent Rq{z) = i(z — Lq) 1 . Here, we 
have used the identities 



Tr L (c /fc p L ) = Tv L {c\ k p L ) 



0. 



P T L V P T = P T IiP T = 0, 
[P,L ] = [Q,L }=0, 
UP = PL = 0. 



(4.12) 

(4.13) 
(4.14) 
(4.15) 



Equation ((0|) follows from Eq. (fl|), while Eq. ( |4.14| ) 
holds because Hq neither mixes the QDS with the leads 
nor does it change the diagonal elem ents or the particle 
num ber o f a state. Finally, Eq. (4.15) can be shown with 
Eq. ( 4.14 ) and using that P contains Pjj. For an expan- 
sionin the small perturbation Ly in Eqs. (4.7), (4.£) and 
(4.11) we use the von Neumann series 



z — QLQ z — Lo — QLyQ 

oo 

= -iR (z)Q [-iL v Ro(z)Q] n . (4.16) 



n=0 



B. Master Equation 



Using Eqs. (4.3), (14), and (|4.10D the diagonal part of 
the reduced density matrix ps(z) — PdPn^lp(z) can 
now be written as 



p s (z) = Tr L PR(z)Pp = 



z-S(z) 



Ps- 



(4.17) 



This equation leads to ps(z) = —izps(z) — ps = 
—iYi(z)ps{z). The probability p n {z) = (n\ps(z)\n) for 
the QDS being in state |n) then obeys the equation 



Pn( z ) = X! W nm(z)p m (z), 



(4.18) 



W nm (z) = -iTr s p„^(z)p m = -iT, nn \ mm (z), (4.19) 

with p n = \n){n\, which is a closed equation for the 
density matrix in the subspace defined by P (with fixed 



N) . In the cotunneling regimecll, the sequ ential tunneling 
contribution (second order in Ly) to Eq. (4.19 ) vanish es. 
The leading contribution [using Eqs. (4.11) and ( 4.16| )] is 
of fourth order in Ly, 



W nm = Tr p n (L v QR ) 3 L v p 

mPh ■ 



(4.20) 



Note that since we study the regime of small frequen- 
cies Rez = uj <C H^oQH ~ l-^n — E m \, where m ^ n, 
we can take the limit u — > here. In addition to this, 
we have assumed fast relaxation in the leads and have 
taken the Markovian limit z — irj — > 0, i.e. we have re- 
placed Wnmj z) in Eq. fl4.19|) by W nm = lim 2 _ W nm {z) 



in Eq. ( 4.20| ) . The trac e of p is preserved under the 
time evolution Eq. ( 4. lg ) since J2 n W n m has the form 
TrP N L v A = Tr [V, A] - Tr Q N [V, A] where the first term 
vanishes exactly and the second term invloving Qn = 
1 — Pn is O(k). After some calculation, we find that 
W nm is of the form 



Wr, 



(4.21) 



with w nm > for a ll n and to. Substituting this equa- 
tion into Eq. ( 4. IS ) we can rewrite the master equa- 
tion in the manifestly trace-preserving form p n {z) — 

(z)}, or in real time, 



Pn{t) = ^2 [ W nraPrn{t) 



iPn(t)} 



(4.22) 



This "classical" master equation describes the dynamics 
of the QDS, i.e. it describes the rates with which the 
probabilities p n for the QDS beinjg_in state \n) change. 
After some algebra (retaining onlycll O(k ), cf. App. |l|), 
we find 



Wnm — w nm + w nm + w nmi 

where (in the cotunneling regime) 



W nm w rwn(2, 1), W nm — w nm{^-j 2), 

w nm = W nm(lJ), 
1=1,2 



(4.23) 

(4.24) 
(4.25) 



with the "golden rule" rate from lead I to lead I', 
w nm (l',l) = 2tt^ |(n|(£>J, A')|m)| 2 

x5(E m -E n -Ap w )p L ,m- (4.26) 

In this expression, Apw = pi — py denotes the chemical 
potential drop between lead I and lead I', and pL,m = 
(to | /?l | to.). We have defined the second order hopping 
operator 

(D},D V ) = DjR Q D v + d v r q d\ 

= -{D\D v +D v D\), (4.27) 
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where Di is given in Eq. (Q, Di = f"^ Dt(t)dt. 

Note, that (Dj, Dp) is the amplitude of cotunneling from 
the lead I to the lead V (in particular, we can write 
B = -(Dl,D2), see Eq. (|3.15|)). The combined index 
m = (m, fh) contains both the QDS index m and the 
lead index fh. Correspondingly, the basis states used 
above are |m) = \m)\m) with energy E m — E m + E m , 
where \m) is an eigenstate of H$ + -ffij nt with energy E m , 
and |m) is an eigenstate of Hl + Hf^ — ^2nPiNi with 
energy E m . The terms w^ m account for the change of 
state in the QDS due to a current going from lead 1 to 
2 (2 to 1). In contrast to this, the cotunneling rate w" m 
involves either lead 1 or lead 2 and, thus, it does not con- 
tribute directly to transport. However, w^ m contributes 
to thermal equilibration of the QDS via particl e-hole ex- 
citations in the leads and/or QDS (see Sees. VIA and 
VTB|). 



C. Stationary State 

In order to make use of the standard Laplace trans- 
form for finding the stationary state p of the system, we 
shift the initial state to to = and define the stationary 
state as p = lim^oo p(t) = lim^oo e~ zLt po. This can be 
expressed in terms of the resolvent, 



p = —i lim zR{z)po, 



(4.28) 



using the property limt—.oo f(t) = —i lim z ^o z f( z ) °f the 
Laplace transform. The stationary state ps of t he QDS 
can be obtained in the same way from Eq. (4.17), 



Ps 



lim — 

z^O z - 



E(z) 



Ps- 



(4.29) 



Multiplying both sides with z — E(z) and taking the limit 
z — » 0, we obtain the condition 



EoPs = 0, 



(4.30) 



where Eo = liniz->o E(z). Using Eq. ( 4.19] ), this condition 
for the stationary state can also be expressed in terms of 
W nm , 



WmnPn) = 0, (4.31) 



which is obviously the s tationarity condition for the mas- 
ter equation, Eq. ( 4.22 ). 



D. Average Current 



The expectation value Ii(t) — Tr Iip(t) of the current 



in lead I [Eq. (2.7)] can be obtained via its Laplace 
transform 

Il(z) = Tr i lP (z) = Tr I l (P + Q)R(z)Pp , (4.32) 



wher e we have inserted P+Q = 1 and used Eqs. (4.3) and 
(4.4) for p(z). According to Eq. ( 4.1 3| ) the first term van - 
ishes. The second term can be rewritten using Eqs. (4.7) 
and (4.17), with the result 



k{z) = Tr/ ; Q 



-QL v p s (z)p L 



z-QLQ 

= Tv S W I {z)ps{z) =J2 W nm(^Pm(z). (4.33) 
nm 

Using the projector method, we have thus managed to 
express the expectation value of the current (acting on 
both the QDS and the leads) in terms of the linear su- 
peroperator W 1 which acts on the reduced QDS density 
matrix p$ only. Taking z — > in Eq. (4.33), the average 



current in lead / in the stationary limit becomes 



= lim Tt IiQ — 

z^0 Z - 



1 



QLQ 



QLvPsPl- 



(4.34) 



Up to now this is exact, b ut ne xt we use again the 
pertur batio n expansion Eq. ( 4.16| ). In the cotunneling 
regimeE3c3, i.e. away from resonances, the second-order 
tunneling current 



l\ 2) = -iTTliR L v psPL 



(4.35) 



is negligible [O(k)], and the leading contribution is the 
cotunneling current 



r( 4 ) 



iTr It(QR L v ) 3 p sPL . 



(4.36) 



After further calculation we find in leading order (cf. 
Ap P . § 



(4.37) 
(4.38) 



where w^ m a re de fined in Eq. (4.24). Note again that 
w^m m Eq. ( 4.25| ) does not contribute to the c urren t 
directly, but indirectly via the master equation Eq. (4.31) 
which determines p m (note that p m is non-perturbative 
in V). We finally remark that for Eqs. (^3^)-( |4.37| ) we 
do not invoke the Markovian approximation. 



E. Current Correlators 

Now we study the current correlators in the stationary 
limit. We let to — > — oo, therefore p(t = ) — > p. The 
symmetrized current correlator [cf. Eq. (|2.7j)], 



S w (t) ^ReTr SIi(t)SI v p, 



(4.39) 



where SIi = Ii — Ii, can be rewritten using the cyclic 
property of the trace as 



S w (t) =ReTr5I l e- ltL 5I l ,p 



(4.40) 
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where e~ ltL acts on everything t o its right. Taking the 
Laplace transform and using Eq. ( 4.28 ) for the stationary 
state p, we obtain 

S w {z)= Km Re(-iz')Tr5IiR(z)5Ii-R(z')Pp , (4.41) 

Q = l 



where z — uj + ir\ and r\ — > 0+. We insert P 
twice and use Eq. ( 4 . 1 2j ) with the result 



S w (z) = S&[z) + &-(i/z)I l I v 



(4.42) 



where S% = S^® + S% P '. We further evaluate the con- 



tributions to Sw(z) using Eqs. (4.7) and (4.29), and we 
obtain 



where Rq ■■ 



qQQ 



■- ReTr IiRqL v PR(z)PIi>RqLvp, 
lim z ^o(z — QLQ) -1 , and 

= —KeTr IiRqLvQRoIi' RoLyp 
—ReTi IiRoIi'QRoLvRoLvp, 
■ — KeTr IiRqLvQRqLvRqIi'P- 



(4.43) 



(4.44) 
(4.45) 



While Sft, (z) as given in Eq. (| 4-43| ) is a non-perturbative 
result, we have used Eq. (4.16) to find the leading con- 



tribution in the tunneling amplitude Tik p for Sffi and 



S% P in Eqs. (03) and (OJ). Also note that QR(z)Q 



was replaced by QR Q in Eqs. (|L4J) and Q4 4E 
w < l-En - -Em | for n ^ m and therefore ,S$ y and S$ P 
do not depend on z, i.e. they do not depend on the fre- 
quency uj. 

In order to analyze Eq. ( 4.43] ) further, we insert the 
resolution of unity ^2 m p m = Is next to the P operators 
in Eq. flQjjj ) with the result S£ = Sg 2 = -Sf 2 = -S£ 
where 



Sf 1 = AS +(»/*)# 
with the non-Poissonian part 



(4.46) 



AS(z) = e 2 ^ w I nm 5p mnl {z)w I , m ,p ml . (4.47) 



The conditional density matrix is defined as 

<fyWi(*) = p nm (z) - (i/z)p n , 
p nm {z) = Tr p n R(z)p m p L . 



(4.48) 
(4.49) 



Eq. (4.17) shows that p nm (z ) must be a solution of 
the master equation Eq. (4.22) for the initial condition 
As(0) = Pm, or p„(0) = 8 nm . We now turn to the remain- 



We find = = -5? 2 = = S Q , where 

S Q = e 2 ^K| m + w- m )p m . (4.51) 



Fina lly, w e can combine Eqs. ( |4.47 ) and (4.51), using 
Eq. (]4.42|) and we obtain the final result for the current 
correlators, 

Sn(w) = S 22 (oj) = -S 12 (uj) = -S 21 (oj) = S(oj), (4.52) 
S(u) = e 2 ^2(w+ m + w~ m )p m + AS(w), (4.53) 



AS(w) = e 2 ^ 1»nm S Pmn> {^W n , m ,p m - , 



(4.54) 



n,m,n' ,m' 



where 8p nm (w) = p„ m (w) - 27r8(uj)p n . Here, p nm (oj) 
is the Fourier-transformed conditional density matrix, 
which is obtained from the symmetri zed so lution p n (t) = 
p n (—t) of the master equation Eq. ( 4.22 ) with the ini- 
tial condition p n (0) = S nm . Not e tha t p nm (uj) is re- 



lated to the Laplace transform Eq. ( 4.49 ) via the relation 

For a few-level QDS, SE ~ Ec, with inelastic cotun- 
neling the noise will be non-Poissonian, since the QDS 
is switching between states with different currents. An 
explicit result for the noise in this case can be obtained 
by making further assumptions about the QD S and the 
coupling to the leads, and then evaluating Eq. ( 4.54 ), see 
the following sections. For the general case, we only es- 
timate AS". The current is of the order / ~ ew, with w 
some typical value of the cotunneling rate w nm , and thus 
51 ~ ew. The time between switching from one dot-state 
to another due to cotunneling is approximately r ~ it; -1 . 
The correction AS to the Poissonian noise can be esti- 



mated as AS ~ SI 2 t 



which is of the same order 



ing contribution Sjj, to Sw(z) in Eq. ( 4.42 ). The Fourier u } ar H- t = H 



as the Poissonian contribution el ~ e 2 w. Thus the cor- 
rection to the Fano factor is of order unity. In contrast to 
this, we find that for elastic cotunneling the off-diagonal 
rates vanish, w nm oc 5 nrn , and therefore Sp nn = and 
AS = 0. Moreover, at zero temperature, either w^ n or 
w~ n must be zero (depending on the sign of the bias Ap). 
As a consequence, for elastic cotunneling we find Poisso- 
nian noise, F = S(0)/e\I\ = 1. 

In sum mary, we have derived the mas ter equation, 
Eq. ( [4.22| ), the stationar y stat e Eq. ( |4.29| ) of the QDS, 
the averag e curr en t, Eq . (4.37), and the current correla- 
tors, Eqs. fl4.52| )- fl4.54p for the QDS system coupled to 
leads in the cotunneling regime under the following as- 
sumptions. (1) Strong cotunneling regime, w- m -c w, i.e. 
negligible intrinsic relaxation in the QDS compared to 
the cotunneling rate; (2) the weak perturbation V is the 
only coupling between the QDS and the leads, in partic- 



transform Sjj7(u)) of the noise spectrum can be obtained 
from its Laplace transform SjfF(z) by symmetrizing the 
latter, 



Sf l T(oj) = S^(uj) + Sfif(-oj). 



LT / 



(4.50) 



s -r Hf\ where Hp acts on the QDS 
and H™ 1 on the leads only; (3) no quantum correlations 
(neither between the QDS and the leads nor within the 
QDS or the leads) in the initial state, po — Ppo; (4) no 
degeneracy in the QDS, E n =^ E m for n ^ m; (5) small 
frequencies, to -C \E m — E n \. For the master equation 
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Eq. ( 4.22| ) (but not for the other results) we have ad- 
ditionally used the Markovian approximation, assuming 
fast relaxation in the leads compared to the tunneling 
rate. 



V. COTUNNELING THROUGH NEARLY 
DEGENERATE STATES 

Suppose the QDS has nearly degenerate states with 
energies E n , and level spacing A nm = E n — E rn , which is 
much smaller than the average level spacing SE. In the 
regime, Ap, k B T, A nm <C SE, the only allowed cotunnel- 
ing processes are the transitions between nearly degener- 
ate st ates. The noise power is given by Eqs. (4.53) and 
(4.54), and below we calculate the correlation correction 
to the noise, AS. To proceed with our calculation we 



rewrite Eq. ( |4.22| ) for Sp(t) (see Eq. ( |4.48| )) as a second- 
order differential equation in matrix form 



5p(t) = W 2 Sp(t), Sp(0) = l-p, 



(5.1) 



where W is defined in Eq. (4.21). We solve this equation 
by Fourier transformation, 



Sp{uj) 



2W 



W 2 



L0 2 1' 



(5.2) 



where we have use d W p — 0. We substitute 5p from this 
equation into Eq. ( 4.54 ) and write the result in a compact 
matrix form, 



AS(w) 



2W 



w 



W 2 +u 2 l 



w 1 p 



(5.3) 



This equation gives the formal solution of the noise prob- 
lem for nearly degenerate states. As an example we con- 
sider a two-level system. 

Using the detailed balance equation, W2\P\ = 
W\2p2, we obtain for the stationary probabilities p\ — 
w\il {w\i + w 2 i), and p 2 = W2i/(u , i 2 + w 2 i)- From 
Eq. ( |4.37D we get 



Wl2{w[\ + U>2i) + W2l(W22 + w{ 2 ) 



Wl2 + W21 



(5.4) 



A straightforward calculation with the help of Eq. (5.2) 
gives for the correction to the Poissonian noise 



. , 2e 2 (wL + wL — wL — wh) 

(Wi2 + W2l)[UJ z + (W12 + W2l) J 

(1 - 2)] . 



X [wl 1 Wi2W21 + W^U) 2 ! 



(5.5) 



In particular, the zero frequency noise AS(0) diverges if 
the "off-diagonal" rates w nm vanish. This divergence has 
to be cut at uj, or at the relaxation rate u>i n due to cou- 
pling to the bath (since u>i2 in this case has to be replaced 
with u>i2 + Win)- The physical origin of the divergence 
is rather transparent: If the off-diagonal rates u>i2,u>2i 



are small, the QDS goes into an unstable state where it 
switches between states 1 and 2 with different currents 
in general.eS The longer the QDS stays in the state 1 or 
2 the larger the zero-frequency noise power is. However, 



if 



21 



u 22 



J l2i 



then AS(u) is suppressed to 



0. For instance, for the QDS in the spin-degenerate state 
with an odd number of electrons AS(u>) = 0, since the 
two states | f) and | |) are physically equivalent. The 
other example of such a suppression of the correlation 
correction AS to noise is given by a multi-level QDS, 
SE <C Ec, where the off-diagonal rates are small com- 
pared to the diagonal (elastic) rates.Ea Indeed, since the 
main contribution to the elastic rates comes from tran- 
sitions through many virtual states, which do not par- 
ticipate in inelastic cotunneling, they do not depend on 
the initial conditions, w\i = u>2 2 , and cancel in the nu- 



merator of Eq. (5.5), while they are still present in the 
current. Thus the correction AS/I vanishes in this case. 
Further below in this section we consider a few-level QDS, 
5E ~ Ec, where AS ^ 0. 

To simplify further analysis we consider for a moment 
the case, where the singularity in the noise is most pro- 
nounced, namely, oj = and |Ai 2 | <C A/i, fcgT, so that 
= and wi 2 = W2i- Then, from Eqs. (5.4) and 
.50 we obtain 



1 = \{h + h) , 4 = e^ wlnn , 



m=l,2 



A5(0) 



4wi2 



(5.6) 



(5.7) 



where I n is the current through the n-th level of the 
QDS. Thus in case |Ai 2 | < Ap,,k B T the following 
regimes have to be distinguished: (1) If k B T < A/z, then 
I n oc Ap, w\2 oc Ap, and thus both, the total current 
I = e~ 1 G dA(i, and the total noise S = FGoAp are 
linear in the bias Ap, (here Go is the conductance of 
the QDS). The total shot noise in this regime is super- 
Poissonian with the Fano factor F ~ I/{ewvi) ^1- (2) 
I n th e regime Ap < k B T < F 1 / 2 Ap the noise correction 
(|5.7|) arises because of the thermal switching the QDS 
between two states n — 1,2, where the currents are lin- 
ear in the bias, /„ ~ GoAp/e. The rate of switching 
is w 12 oc k B T, and thus AS - FG D Ap 2 / {k B T). Since 
k B T I Ap < F 1 / 2 , the noise correction AS is the domi- 
nant contribution to the noise, and thus the total noise«£ 
can be interpreted as being a thermal telegraph noise.cS 
(3) Finally, in th e reg ime F 1 l 2 Ap < k B T the first term 
on the rhs of Eq. (4.53) is the dominant contribution, and 
the total noise becomes an equilibrium Nyquist noise, 
S = 2Gok B T. 

We notice that for the noise power to be divergent the 
off-diagonal rates w\2 and u>2i have to vanish simultane- 
ously. However, the matrix w nm is not symmetric since 
the off-diagonal rates depend on the bias in a different 
way. On the other hand, both rates contain the same ma- 
trix e lemen t of t he co tunneling amplitude (Dj, Dy), see 
Eqs. ( 4.26| ) and (4.27). Although in general this matrix 



11 




FIG. 2. Double-dot (DD) system containing two electrons 
and being weakly coupled to metallic leads 1, 2, each of which 
is at the chemical potential /ii , fi2 ■ The tunneling amplitudes 
between dots and leads are denoted by T. The tunneling (td) 
between the dots results in a singlet-triplet splitting J ~ t d /U 
with the singlet being a ground stateO The tunneling path 
between dots and leads 1 and 2 forms a closed loop (shown 
by arrows) so that the Aharonov-Bohm phase <f> will be accu- 
mulated by an electron traversing the DD. 



element is not small, it can vanish because of different 
symmetries of the two states. To illustrate this effect we 
consider the transport through a double-dot (DD) sys- 
tem (see Ref. ^l] for details) as an example. Two leads 
are equally coupled to two dots in such a way that a 
closed loop is formed, and the dots are also connected, 
see Fig. ^. Thus, in a magnetic field the tunneling is 
described by the Hamiltonian Eq. ( |2.3| ) with 

D l = Y,Ti j 4 s djs, 1,3=1,2, (5.8) 

T n = T 22 = T* 2 = T" 2 *i = e I0/4 T , (5.9) 

where the last equation expresses the equal coupling of 
dots and leads and (f> is the Aharonov-Bohm phase. Each 
dot contains one electron, and weak tunneling td between 
the dots causes the exchange splitting^ J ~ t d /U (with 
U being the on-site repulsion) between one spin singlet 
and three triplets 



1 



V2 



■\dLdf 



\S) 

|T + ) = 4 T 4 T |0), \T_) = d\ i d\ i \Q) 



(5.10) 



In the case of zero magnetic field, = 0, the tunnel- 
ing Hamiltonian V is symmetric with respect to the ex- 
change of electrons, 1^2. Thus the matrix element of 
the cotunneling transition between the singlet and three 
triplets (S\V(E~ Ho)^ 1 ]^^), i — 0,±, vanishes because 
these states have different orbital symmetries. A weak 
magnetic field breaks the symmetry, contributes to the 
off-diagonal rates, and thereby reduces noise. 

The fact that in the perturbation V all spin indices are 
traced out helps us to map the four-level system to only 
two states \S) and |T) classified according to the orbital 
symmetry (since all triplets are antisymmetric in orbital 



Q.=0.l 




0.4 0.6 0.8 1.0 

<j)/27I 



FIG. 3. The Fano fact or F = S(l j)/L with the noise 
power S(ui) given in Eqs. (4.53) and (5.12), and wit h th e 
current through the DD, I, given in Eqs. (5.4) and (5.11), 



is plotted as a function of the Aharonov-Bohm phase (j> for 
the normalized bias v = A/i/J = 2 and for four different 
normalized frequencies SI = u/[G(2A/j, — J)] = 0.1, 0.25, 0.5, 
and 1. Inset: the same, but with fixed frequency £1 = 0.1, 
where the bias v takes the values 1.5, 3, and oo. 



space) . In Appendix |c] we derive the mapping to a two- 
level system and calculate the transition rates w+ m and 
w nm ( n , m ~ 1 f° r a si nglet and n, m = 2 for all triplets) 
using E qs. ( 4.26 ) and ( 1.27 ) with the operators Di given 
by Eq. (|5.8|). Doing this we obtain the following result 



= 0, 



/T 2 



(1 + cos 4>)A^ 



(1-cos 0)(A/X+J) 



3(1 — cos (f>)(Afx—J) 3(l+cos0)A/x 



, (5.11) 



which holds close to the sequential tunneling peak, A_ <C 
A + ~ U (but still A_ > J, An), and for Afj, > J. We 
substitute this equation into the Eq. ( |5.5| ) and write the 
correction AS(u>) to the Poissonian noise as a function 
of normalized bias v = Aji/ J and normalized frequency 
Q = ecj/[G(2Afi — J)] 



AS(lj) = 6eGJ 



{v 2 - 1)[1 + (v- l)cos0] 2 (l 
(2v- l) 3 [n 2 + (1 -cos<; 



COS ( 



(5.12) 



where G = 7re(i/T 2 /A_) 2 is the-iconductanc e of a single 
dot in the cotunneling regime.Ell From Eq. (5.12) it fol- 
lows that the noise power has singularities as a function 
of uj for zero magnetic field, and it has singularities at 
<\> = 2mn (where m is integer) as a function of the mag- 
netic field (see Fig. |^). We would like to emphasize that 
the noise is singular even if the exchange between the 
dots is weak, J <C A/i. Note however, that our classical 
approach, which neglects the off-diagonal elements of the 
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density matrix p(t), can only be applied for weak enough 
tunneling, w nm <C J ■ In the case A/i < J the transition 
from the singlet to the triplet is forbidden by conserva- 
tion of e nergy, = 0, and we immediately obtain from 
Eq. (5.5) that AS(uj) = 0, i.e. the total noise is Poisso- 
nian (as it is always the case for elastic cotunneling). In 
the case of large bias, A/z ^> J, two dots contribute in- 
dependently to the current I — 2e~ 1 GA/z, and from Eq. 
( |5.12 ) we obtain the Fano factor 



F 



3 cos<Wl 



COS( 



Q 2 + (1 - cost 



A/i > J. 



(5.13) 



This Fano factor controls the transition to the telegraph 
noise and then to the equilibrium noise at high tempera- 
ture, as described above. We notice that if the coupling 
of the dots to the leads is not equal, then w° m ^ serves 
as a cut-off of the singularity in AS(uj). 

Finally, we remark that the Fano factor is a periodic 
function of the phase <j> (see Fig. ||); this is nothing 
but an Aharonov-Bohm effect in the noise of the co- 
tunneling transport through the DD. However, in con- 
trast to the Aharonov-Bohm effect in the cotunneling 
current through the DD which has been discussed earlier 
in Ref. |l|, the noise effect does not allow us to probe 
the ground state of the DD, since the DD is already in a 
mixture of the singlet and three triplet states. 



VI. COTUNNELING THROUGH CONTINUUM 
OF SINGLE-ELECTRON STATES 

We consider now the transport through a multi- 
level QDS with SE <C Ec- In the low bias regime, 
Aa <C (ciE.Ec) 1 ^ 2 , the elastic cotunneling dominates 



transport,E£l and according to the results of Sees. |IV| and 
^ the noise is Poissonian. Here we consider the oppo- 
site regime of inelastic cotunneling, A/i ^ (5E Ec) 1 ^ 2 ■ 
Since a large number M of levels participate in transport, 
we can neglect the correlations which we have studied 
in the previous section, since they become a 1/M-effect. 
Instead, we concentrate on the heating effect, which is 
not relevant for the 2-level system considered before. 
The condition for strong cotunneling has to be rewrit- 
ten in a single-particle form, T- m 3> t c , where T- m is the 
single-particle energy relaxation time on the QDS due 
to the coupling to the environment, and r c is the time 
of the cotunneling transition, which can be estimated as 
t c ~ ev^A^/I (where i>u is the density of QDS states). 
Since the energy relaxation rate on the QDS is small, the 
multiple cotunneling transitions can cause high energy 
excitations on the dot, and this leads to a nonvanishing 
backward tunneling, w~ m =/= 0. In the a bsenc e of cor rela- 
t ions between cotunneling events, Eqs. (4.37), ( 4.3§| ) and 
(4.53) can be rewritten in terms of forward and backward 
tunneling currents i+ and 



I = U-I- 



S = e(I++I.), 



(6.1) 



I± 



e 7 , w nmPm , 



(6.2) 



where the transition rates are given by Eqs. (4.24) and 

dli- 
it is convenient to rewrite the currents I± in a single- 
particle b asis . To do so we substitute the rates Eq. ( 4. 26] ) 
into Eq. ( |6.2| ) and n egle ct the dependence of the tunnel- 
ing amplitudes Eq. (2.3) on the quantum numbers k and 
p, Tikp = Ti, which is a reasonable assumption for QDS 
with a large number of electrons. Then we define the 
distribution function on the QDS as 



^^e-e^Trpdtdp 



(6.3) 



and replace the summation over p with an integration 
over e. Doing this we obtain the following expressions 
for T = 



I±=C ± 



C± = 



G1G2 
2ne 3 



1 

A~ 



1 

aZ 



(Am) 3 , 



(6.4) 



Afi 3 



dede'0(e- e'± A/*)/(e)[l - f(e% (6.5) 



where Gi ; 2 = Tre 2 uuu\Ti,2\ 2 are the tunneling conduc- 
tances of the two barriers, and where we have intro- 
duced the function 0(e) = ed(e) with 6(e) being the 
step- function. In particular, using the property Q(e + 
A/i) — 8(e — A/i) — e + A/i and fixing 



de[f(e) 



-e)] = 0, 



(6.6) 



(since I± given by Eq. (6.4) and Eq. ( |6.5| ) do not depend 
on the shift e — > e + const) we arrive at the following 
general expression for the cotunneling current 



/ = A 



G1G2 



1 

AT 



1 



A = 

T 



127re 3 \A + ' A 
1 + 12T/(A/i) 2 , 

dee[f(e) - 6(-e)} > 



(Am) 3 , 



(6.7) 
(6.8) 
(6.9) 



where the value vj^T has the physical meaning of the en- 
ergy acquired by the QDS due to the cotunneling current 
through it. 

We have deliberately introduced the functions C± in 
the Eq. ( |6.4[ ) to emphasize the fact that if the distribu- 
tion /(e) scales with the bias A/i (i.e. / is a function 
of e/A/t), then C± become dimensionless universal num- 
bers. Thus both, the prefactor A [given by Eq. fl6.8|) ] in 
the cotunneling current, and the Fano factor F = S/ (el), 
where S = el + AS h , 



F = 



C+ + C- 
C+-C- 



(6.10) 



take their universal values, which do not depend on the 
bias A/x. We consider now such universal regimes. The 
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first example is the case of weak cotunneling, r; n <C r c , 
when the QDS is in its ground state, /(e) — 0(—e), and 
the thermal ener gy o f the QDS vanishes, T = 0. Then 
A = 1, and Eq. (5.7) reproduces the results of Ref. |2^. 
As we have already mentioned, the backward current van- 
ishes, I- = 0, and the Fano factor acquires its full Poisso- 
nian value F = 1, i n agre ement with our noncquilibrium 
FDT proven in Sec. Ill B . In the limit of strong cotunnel- 
ing, r; n ^> r c , the energy relaxation on the QDS can be 
neglected. Depending on the electron-electron scattering 
time T ee two cases have to be distinguished: The regime 
of cold electrons r ee ^ r c and regime of hot electrons 
Tee <C t c on the QDS. Below we discuss both regimes in 
detail and demonstrate their universality. 



A. Cold electrons 

In this regime the electron-electron scattering on the 
QDS can be neglected and the distrib ution /(e) has to be 
found from the master equation Eq. (4.22). We multiply 
this equation by ^2 p S(e — e p )(n\dj ) dp\n), sum over n 
and use the tunneling rates from Eq. (4.26). Doing this 



we obtain the standard stationary kinetic equation which 
can be written in the following form 



J de'a(e' 



e)f(e')[l-f(e)] 



— J de'a(e 



e')f(e)[l-f{e% (6.11) 



<j{e) = 2A6(e) + ^ 6(e ± A/x), 



(6.12) 



where A = {G\ + 
libration rate iu£. 



G 2 2 )/(2G 1 Go 



> 1 arises from the equi- 
(We assume that if 



„.„. see Eq. fl4.25| ) 
the limits of the integration over energy e are not speci- 
fied, then the integral goes from — oo to +oo.) From the 
form of this equation we immediately conclude that its 
solution is a function of e/A/x, and thus the cold electron 
regime is universal as defined in the previous section. It 
is easy to check that the detailed balance does not hold, 
and in addition er(e) ^= a(—s ). Thus we face a difficult 
problem of solving Eq. ( |6.11 ) in its full nonlinear form. 
Fortunately, there is a way to avoid this problem and to 
reduce the equation to a linear form which we show next . 

We group all nonlinear terms on the rhs of Eq. ( 3.11 ) : 
/ de'a(e' - e)f(e') = h(e)f(e), where h{e) = 
J de' {a(e' ~ e)/(e') + a(e - e')[l - f(e')}}. The trick is 
to rewrite the function h(e) in terms of known functions. 
For doing this we split the integral in h(e) into two inte- 
grals over e' > and e' < 0, and then use Eq. ( |6.6[ ) and 
the property of the kernel a(e) — er(— e) = 2(1 + A)e to 
regroup terms in such a way that h(e) do es n ot contain 
/(e) explicitly. Taking into account Eq. ( |6.9| ) we arrive 
at the following linear integral equation 



J de'a{e' 



2.2 
2.1 
2.0 
< 1.9 
1.8 
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FIG. 4. The prefactor A in the expression (6/7) for 
the cotunneling current characterizes a universal cotunnel- 
ing transport in the regime of weak cotunneling, r m <C t c , 
(A = 1, see Ref. and in the regime of strong cotunnel- 
ing, Ti n 2> t c (A > 1). Here A is plotted as a function of 
G1/G2 (same as a function of G2/G1) for the strong cotun- 
neling, for the cold-electron case, r ee 3> t c (solid line) and 
for the hot-electron case, r ee t c (dotted line). Gi,2 are the 
tunneling conductances of a junctions connecting leads f and 
2 with the QDS. 



[(l + A)(e 2 + 2T) + (A^) 2 ]/(e), 



.13) 



where the para meter T is the only signature of the non- 
linearity of Eq. (|6.ff ) . 

Since Eq. (3.13) represents an eigenvalue problem for 
a linear operator, it can in general have more than one 
solution. Here we demonstrate that there is only one 
physical solution, which satisfies the conditions 

0</(e)<l, /(-oo) = l, /(+oo)=0. (6.14) 

Indeed, using a standard procedure o ne ca n show that 
two solutions of the integral equation (|6.13 ), /1 and /2, 
corresponding to different parameters Ti 7^ T2 should 
be orthogonal, J d e/i(e )/2(— e) = 0. This contradicts 
the conditions Eq. ( |6. 14 ) . The solution is also unique for 
the same T, i.e. it is n ot de generate (for a proof, see A p- 
pendix ^) . From Eq. ( |3.11 ) and conditions Eq. ( |5.14 ) it 
follows that if /(e) is a solution then 1 — /(— e) also sat- 
isfies Eqs. (6.11 ) and (6.14). Since the solution is unique, 
it has to have the symm etry /(e) = 1 — /(— e). 

We sol ve Eq s. (6.13 ) and ( 6.14 ) numerically and use 
Eqs. (6.5) and (3.10) to find that the Fano factor is very 
close to 1 (it does not e xcee d the value F « 1.006). Next 
we use Eqs. ( |6.8[ ) and ( |6.9| ) to calculate the prefactor A 
and plot the result as a function of the ratio of tunnel- 
ing conductances, G1/G2, (Fig. 4, solid line). For equal 
coupling to the leads, G\ = G2 , the prefactor A takes its 
maximum value 2.173, and thus the cotunneling current 
is approximately twice as large compared to its value for 
the case of weak cotunneling, r in <C t c . A slowly de- 
creases with increasing asymmetry of coupling and tends 



14 



to its minimum value A = 1 for the strongly asymmetric 
coupling case GxjGzorGijGx S> 1. 



B. Hot electrons 

In the regime of hot electrons, r ee <C r c , the dis- 
tribution is given by the equilibrium Fermi function 
Jf{e) = [1 + exp(e/fcBT'e)]~ 1 , while the electron temper- 
ature T e has to be found self-consistently from the kinetic 
equation. Eq. (6.11) has to be modified to take into ac- 
count electron-electron interactions. This can be done 
by addin g the electron collision integral I e e(e) to the rhs. 
of (6.11). Since the form of the distribution is known 
we need only the energy balance equation, whic h can be 
derived by multiplying the modified equation ( |6.11 ) by 
e and integrating it over e. The contribution from the 
collision integral I ee (e) vanishes, because the electron- 
electron scattering conserves the energy of the system. 
Using the symmetry /f(e) = 1 — If{—£) we arrive at the 
following equation 



j J dede'f F (e')[l 



f F {e)]a{e' -e)e = 0. (6.15) 



Next we regroup the terms in this equation such that it 
contains only integrals of the form defp{s){- ■ ■)■ This 
allows us to get rid of nonlinear terms, and we arrive at 
the following equation, 



Jdee 3 [f F (e) 



6(-e)] 



3T 2 



(Am) 4 

3(1 + A) 



(6.16) 



which holds also for the regime of cold elect rons. Fi- 
nally, we calculate the integral in Eq. ( |6.16 ) and ex- 
press the result in terms of the dimcnsionless parameter 
a = A^i/k B T e , 



a = 7T [8(1 + A)/5] 



1/4 



(6.17) 



Thus, since the distribution again depends on the ratio 
e/A/i, the hot electron regime is also universal. 

The next step is to substitute the Fermi di strib ution 
function with the temperature given by Eq. ( |6.17| ) into 
Eq. ( |6.5[ ). We calculate the integrals and arrive at the 
closed analytical expressions for the values of interest, 



(6.18) 
(6.19) 




2(1 + A) 



F = 1 



12 



2-7T 2 + cr 



where again A = (G\ + G|)/2GiG 2 > 1. It turns out 



that similar to the case of cold electrons, Sec. VIA, the 
Fano factor for hot electrons is very close to 1 (namely, 
it does not exceed the value F w 1.007). Therefore, we 
do not expect that the super-Poissonian noise considered 



in this section (i.e. the one which is due to heating of 
a large QDS caused by inelastic cotunneling through it) 
will be easy to observe in experiments. On the other 
hand, the transport-induced heating of a large QDS can 
be observed in the cotunneling cur rent t hrough the pref- 
actor A, which according to Eq. 

/5?4 



takes its maxi- 
mum value A = 1 + ^5/4 w 2.118 for Gi = G 2 and 
slowly reaches its minimum value 1 with increasing (or 
decreasing) the ratio Gi/G 2 (see Fig. 4, dotted line). Sur- 
prisingly, the two curves of A vs Gi / G 2 for the cold- and 
hot-electron regimes lie very close, which means that the 
effect of the electron-electron scattering on the cotunnel- 
ing transport is rather weak. 



VII. CONCLUSIONS 

The physics of the noise of cotunneling is discussed in 
the Introduction. Here we give a short summary of our 
results. 

In Sec. [fill , we have derived t he no n-equ ilibri um FDT, 
i.e. the universal relations Eqs. (3.10) and (3.21) between 
the current and the noise, for single-barrier junctions and 
for QDS in the weak cotunneling regime, respectively. 
Taking the limit T, lo — *■ 0, we show that the noise is 
Poissonian, i.e. F = 1. 
Sec. 



In 



IV, we have derived the master equation, 



Eq. ( |4.22D , the stationar y stat e Eq. ( (4.29| ) of the QDS, 



the average current, Eq. (4.37), and the current correla- 
tors, Eqs. (4.52)-( i~5l| ) for a non-degenerate QDS system 
(E n ^ E m , n ^ m) coupled to leads in the strong cotun- 
neling regime w m « w at small frequencies, to <C A mn . 
In contrast to sequential tunneling, where shot noise is 
either Poissonian (F = 1) or suppressed due to charge 
conservation (F < 1), we find that the noise in the inelas- 
tic cotunneling regime can be super-Poissonian (F > 1), 
with a correction being as large as the Poissonian noise 
itself. In the regime of elastic cotunneling F — 1. 

While the amount of super-Poissonian noise is merely 
estimated at the end of Sec. IV, the noise of the cotun- 



neling current is calculated for the special case of a QDS 
with nearly degenerate states, i.e. A nm -C SE, in Sec. [v|, 
where we a pply our results from Sec. |Tv| . The general so- 
lution Eq. (5.3) is further analyzed for two nearly degen- 
erate levels, with the result Eq. (5.5). More information 



is gained in the specific case of a DD coupled to lea ds, 
where we determine the correction to noise Eq. ( 5.12| ) as 
a function of frequency, bias, and the Aharonov-Bohm 
phase threading the tunneling loop, finding signatures of 
the Aharonov-Bohm effect in the cotunneling noise. 
Finally, in Sec. VI, another important situation is stud- 



ied in detail, the cotunneling through a QDS with a con- 
tinuous energy spectrum, SE <C A/i <C Ec- Here, the 
correlation between tunneling events plays a minor role 
as a source of super-Poissonian noise, which is now caused 
by heating effects opening the possibility for tunneling 
events in the reverse direction and thus to an enhanced 
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noise power. In Eq. (6.10), we express the Fano factor 
F in the continuum case in term s of the dimensionless 
numbers C±, defined in Eq. (|6 . 5[) , which depend on the 
electronic distribution function f(e) in the QDS (in this 
regime, a description on the sing le-electron level is ap- 
propriate). The curren t Eq . ( |6.7|) is expressed in terms 
of the prefactor A, Eq. (6.8). Both F and A are then cal- 
culated for different regimes. For weak cotunneling, we 
immediately find F = 1, as anticipated earlier, while for 
strong cotunneling we distinguish the two regimes of cold 
(r ee ^> t c ) and hot (r ee <C t c ) electrons. For cold clec- 



function of operators Di an d Dj , then O(oo) = (see 



trons, we derive the linear integral equation Eq. (6.13) for 
f(e) which is shown to have a unique solution, and which 
is solved numerically. We find that the Fano factor is very 
close to one, 1 < F < 1.006, while A is given in Fig. 4. 
For hot electrons, /(e) is the equili brium Fermi distribu- 
tion, and the Fano factor Eq. ( 5.19 ) and A [Eq. ( |6.18] )and 
Fig. 4] can be computed analytically. Again, the Fano 
factor is very close to one, 1 < F < 1.007, which leads us 
to the conclusion that heating will hardly be observed in 
noise, but should be well measurable in the cotunneling 
current. 
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APPENDIX A: 



the discussion in Sec. [II B ). Next, the currents can be 
represented as the difference and the sum of I\ and I 2 , 



h = (i 2 -A)/2 = i e(A t -A)/2, 
I = (Ji + J 2 )/2 = ie(yt - Y)/2 , 



(Al) 
(A2) 



where X = D 2 + D\ , and Y = D x + D 2 . While for the 
perturbation we have 



v = x + x f = y + y f . 



(A3) 



First we concentrate on the derivation of Eq. (3.14) and 
redefine the average current Eq. (2.7) as I = Id (which 
gives the same result anyway, because the average num- 
ber of electrons on the QDS does not change I a — 0). 

To pr oceed with our derivation, we make use of 
Eq. ( p.q) and expand the current up to fourth order in 
Tikp- 

t 

I = i [ dt [ dt'{i d V{t)V{t')V{t')) 



-00 —00 





i / dt(VI d V(t)V(t)) + c.c. 



(A4) 



Next, we use the cyclic property of trace to shift the 
time dependence to Id- Then we complete the integral 
over time t and use Id(oo) — 0. This procedure allows us 
to combine first and second term in Eq. (A4), 



In this Appendix we present the derivation of 
Eqs. ( |H4| ) and fTl6| ). First we would like to mention 
that the operator B in these equations is just the second- 
order tunneling amplitude, which also appears in the 
tunneling Hamiltonian after the Schrieffer- Wolff transfor- 
mation. Therefore, one might think that the Schrieffer- 
Wolff tran sform ation is the most simple way to derive 
Eqs. ( 3.14 ) and ( [3.16; ). On the other hand, it is obvious 
that the Schrieffer- Wolff procedure being a unitary trans- 
formation gives exactly the same amount of terms in the 
fourth-order expression for the current and noise as that 
of the regular perturbation expansion. The Schrieffer- 
Wolff procedure is useful in the Kondo regime where the 
energy scale is given by the Kondo temperature Tk and 
where the i?-terms in the Hamiltonian lead to a diver- 
gence for T < Tk, while the other terms can be treated 
by perturbation theory (see Ref. [52| ). In our cotunneling 
regime such a divergence does not exist (since the QDS 
is weakly coupled to leads, i.e. A/^, fc^T 3> ^bTr-), and 
we have to analyze all contributions. We do this below 
using perturbation theory. 

In order to simplify the intermediate steps, we use the 
notation 0(t) = J_ dt'0{t') for any operator O, and 
O(0) = O. We notice that, if an operator O is a linear 



I = -i dt([I d V + VI d )V{t)V{t)) + , 



(A5) 



Now, using Eqs . (Al) and ( ]A3|) we replace opera- 
tors in Eq. (|A5|) with X and A T in two steps: I = 
e J ^ dt([X^X^-XX]V(t)V(t))+h.c, where some terms 
cancel exactly. Then we work with V{t)V(t) and notice 
that some terms cancel, because they are linear in c;^ 
and cj k . Thus we obtain 



I=e J dt{[X^X^ - XX] 

— OQ 

x [XHt)XHt) + X(t)X(t)}) + c.c. (A6) 

Two terms XXXX and X^X^X^X^ describe tunneling 
of two electrons from the same lead, and therefore they do 
not contribute to the normal current. We then combine 
all other terms to extend the integral to +00, 



/ = e 



dt{X\t)X\t)XX - XXX\t)X\t)) (A7) 
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Finally, we use f^dtXtyXty) = - dtX(t)X(t) 

(since X(oo) = 0) to get Eq. ( 533 ) with B = XX. Here, 

i — i — 

again, we drop terms D[D[ and -D2-D2 responsible for 
tunneling of two electrons from the same lead, and ob- 
tain B as in Eq. (3.15). 



Next, we derive Eq. ( 3.16 ) for the noise power. At 
small frequencies u> <C A± fluctuations of I s are sup- 
pressed because of charge conservation (see below), and 



we can replace I2 in the correlator Eq. (2.7) with 
Id- We expand S(u>) up to fourth order in T)^, use 
J-x dt id(t)e ±lut = 0, and repeat the steps leading to 



Eq. (A5). Doing this we obtain, 



dtcos(u;t)([V{t),i d (t)}[V,i d ]) . (A8) 



Then, we replace V and Id with X and X'. We again 
keep only terms relevant for cotunneling, and in addi- 
tion we neglect terms of order uj/A± (applying same ar- 
gum ents as before, see Eq. (A9)). We then arrive at 



Eq. (3.16) with the operator B given by Eq. (3.15). 
Finally, in order to show that fluctuations of I s are 



suppressed, we replace Id in Eq. (A8) with I s , and then 



use the operators Y and instead of X and X>. In 
contrast to Eq. (A7) terms such as Y^Y^YY do not 
contribute, because they contain integrals of the form 
dt cos(u>t) Di(t)Di'(t) — 0. The only nonzero contri- 
bution can be written as 



S ss (uj) = 



dicos(^)([f t (t),r(t)][r t ,y] 



(A9) 



where we have used integration by part s and the prop- 
erty Y(oo) — 0. Compared to Eq. ( |3. 16 ) this expression 
contains an additional integration over t, and thereby it 
is of order (ui/A±) 2 . 



APPENDIX B: 

We evaluate the matri x elements of the superoperator 
W I {z) given in Eq. ( 1.35 ) w hich are used to calculate the 
average current I; , see Eq. (1.37). The derivation for the 
master equation ( |l.22|) i s very similar. As for the noise, 
the Sf v term Eq. Q4.51p is again obtained in a similar 
way as the current, whereas the S^, t erm Eq. (1.46) is 
different and is analyzed in Sec. [VE. Since W 1 (z) is 



obtained by taking the partial trace over the leads, its 
matrix elements can be expressed as the sum over lead 
indices 



(Bl) 



where n = (n, n), with n and n enumerating the QDS and 
lead eigenstates. For convenience, we will use the eigen- 
states of Ho in this Appendix, and not the eigenstates of 
K as in the main text. Accordingly, here E n = E n + E fl 
are the eigenenergies of Hq. Taki ng th e stationary limit 
z — > 0, using the definition Eq. ( 1.33 ) and introducing 
the projectors p n — |n)(n|, we can write 



W 1 



lim Ti p n IiQ 



1 



z-QLQ 



QL V Pp„ 



(B2) 



Note that while n denotes a free dummy index in 
Eq. (B2), the state |m) is restricted to the subspace 
where PnP-ki 7^ with fixed particle number N on the 
QDS. Expanding this expression in V, we obtain for the 
lowest nonvanishing order (sequential tunneling) the con- 
tribution —i Enm^l^o^vPniJnn to the rate W£ m , which 
can be expressed as 



27re]T(|<n|A|m)| 2 - 



-|(n|D ! V|)| a ) p L ,m<5(A mn ), 



(B3) 



where A mn = E m — E a . Using Eq. (2.3) and assuming 
that Tikp — T is independent of p and k, we obtain the 
expression for the contribution to W 7 { m due to sequential 
tunneling, 

2nvT 2 J2 (\(n\d p \m)\ 2 [1 - fi{A mn )] 
p 

-|(n|4|m)| 2 /i(A„ m )), (B4) 

where // (e) is the Fermi distribution and v the,density of 
states in the leads. In the cotunneling regimeEll, this con- 
tributipa is proportional to k = e~ A / kBT , therefore we 
drop itell and expand W£ m to the next non- vanishing, i.e. 
fourth, order in V. Doing this, we obtain the cotunneling 
contribution 



— i{IiR[)Lv RoQLy RoLvPm)i 



(B5) 



Stepwise evaluation of the operators and superoperators 
in this expression by the insertion of the identity J2i 
leads to 

U^ 1 = {L v R Q L v p m ) Vi 

= ^2 [^ik-Rkj(ivPm)kj - Rik(LvPm)ik.Vkj} , 
k 

(LvPm)ij = V im S m j - V mi S im , (B6) 
where 7y = and similarly for Vy. Note that 

i 

Rw = lim 
= -iP- 



1J v-*) it) - (Ek - Ei) 

' ^{E, - E S ), (B7) 
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where P stands for the principal value. The current is 
obtained from W^ m by multiplying with the full density 
matrix p m and then summing over m and n. By explicit 
evaluation, using the fact that we can choose the basis \n) 
on the QDS such that all expectation values of the form 
(n\d pl dp 2 d p3 d Pi \n} , etc., are real, we find that four out of 
the eight terms in Eq. ( JB6| ) cancel, while the remaining 
four terms contributing to the current can be combined 
into (retaining only O(kP) terms) 

E W nm = -27TlmE [tf (W^ta 
n f 

+ (VR m V) mf (iiRlv) fm \ S(E t - En), (B8) 

where R m = —i¥(Hv — E m )zl. All other (^-function 
contributions vanish in O(k )M3 In the presence of an 
Aharonov-Bohm pha se, when the phases in the tunneling 
amplitudes Eq. (|5.9|) have to be taken into account, we 
again find Eq. ( p38| ) by explicit analysis. We note here 
that exactly the same procedure as above can be applied 
in the derivation of the the master equation and the noise, 
leading to a reduction of terms and finally to the "golden 
rule" expressions Eqs. (4.23) and (4.51). By substituting 
Eqs. ( |2.3| ) and (2.6) for V and Ii, and setting / = 2 for 
concreteness, we finally obtain 

E W nm - 2vr e E [(£>!, D 1 ) mt (D\,D 2 ) tm 



where (D\,D v ) is defined in Eq. ( f4.27|). U sing Eqs. (|4.33| ) 



and (Bl) and the definitions Eqs. (4.24) and (4.26), we 



find for the cotunneling current 



^2=E W n m 

II III 



Pm PL , m 



+ 

nm 



w mn )p m , (BIO) 



which concludes t he de rivation of Eqs. (4.37) and (4.38). 
Note that in Eq. ( 4.26| ) the expression A mn = E m — E n 
is replaced by E m — E n — Apw because there, |n) are 
eigenstates of K (instead of Ho). The current I\ in lead 

1 can be obta ined by interchanging the lead indices 1 and 

2 in Eq. (B9) which obviously leads to I\ = — Ii- 



APPENDIX C: 



In this Appendix we calculate the transition rates 
Eq. ( 4.26| ) for a D P c ouple d to leads with the coupling 



described by Eqs. (5J3) and ( |5.9| ) and show that the four- 
level system in the singlet-triplet basis Eq. (5.10) can 
be mapped to a two-level system. For the moment we 
assume that the indices n and m enumerate the singlet- 
triplet basis, n,m = S, Xo,T+,T_. Close to the sequen- 
tial tunneling peak, A_ <C A + , we keep only terms of 
the form dJRqDi, Calculating the trace over the leads 
explicitly, we obtain at T = 0, 



W nm (l',l) 



2At 



Q(fj,i- fj,i> - A nm ) 



x T*j T\j> Tjtj., Tyj M nm (j, f) , (CI) 

3,3' 

M nm (j,f) - EW d lAil m >H4j' d *j>> , (C2) 



with 0(e) = ed(e), and A nm = 0, ±J, and we have as- 
sumed td -C A_ so that Rq = 1/A_. 

Since the quantum dots are the same we get 
M nm (l,l) = M nm (2,2), and M„ m (l,2) = M nm (2,l). 
We calculate these matrix elements in the singlet-triplet 
basis explicitly, 




(C3) 



(C4) 



{DlD 2 ) mf (DlD 1 ) tm S(A tm ), (B9) follows that 



M(l,l)=o 



M(l,2) 



Assuming now equal coupling of the form Eq. ( |5.9| ) we 
find that for I = I' the matrix elements of the singlet- 
triplet transition vanish (as we have expected, see Sec. 
M). On the other hand the triplets are degenerate, i.e. 
= in the triplet sector. Then from Eq. (CI) it 



= E; 



0. Next, we have 



0(a*2 — fJ-i — A nm ) = 0, since for nearly degenerate states 



we assume Afi > |A„ 
Finally, for w+ = w nrr , 



, and thus w nm - 
(2, 1) we obtain, 



1,21 



J ss 



J ST 



w 



TS 




A/i(l 



(Afj, + J)(l - cos< 



(An- J)(l 



COS( 



0. 

(C5) 
(C6) 
(C7) 



1 + COS (f> 1 + COS <j> 

2 + 2cos0 
2 + 2 cos < 



(C8) 



Next we prove the mapping to a two-level system. 
First we notice that because the matrix w^ T is sym- 
metric, the detailed balance equation for the stationary 
state gives p n /p m = w mn/ w nm = 1' n,m € T. Thus 
we can set p n — * P2/3, for n 6 T. The specific form 
of the transition matrix Eqs. (|C5|-C8) helps us to com- 
plete the mapping by setting (1/3) YL.=2 w im ~* ^12) 



Em=2 



and (1/3) E n,m—2 



SO 



that we get the new transition matrix Eq. (5.11), while 
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the stationary master equation for the new two-level den- 
sity matrix does not change its form. If in addition we set 

(V 3 ) Em=2^Plm(i) -> Sp 12 [t), J2t=2 S Pnl{t) ~> Sp 2 l(t), 

and (1/3) J2n . m=i 5pnm(t) — > 5p22(t), then the master 
equation Eq. ( fl.22| ) for Sp nm (t) and the initial condition 
^Pnm(O) = 5 nm — p n do not change cithe r. Finally, one 



can see that under this mapping Eq. ( 4.54 ) for the correc- 
tion to the noise power AS(w) remains unchanged. Thus 
we have accomplished the mapping of our singlet-triplet 
system to the two-le vel sy stem with the new transition 
matrix given by Eq. (5.11). 



condition Eq. (D2) ensures that the solution of Eq. ( |D5| ) 
is localized, ip(x)\ x ^± oc — and finite everywhere. All 
these requirements can be satisfied only if <p(x) = for all 
x. Indeed, since the function u(x) + (Ap) 2 + 2(1 + A)T 
is positive for all x (we recall that T > 0), then (p is 
a monotonous function, and therefore it cannot be lo- 
calized. In other words, the Schrodinger equation with 
repulsive potential u(x) > does not have localized so- 
lutions. Thus we have pro ven t hat /i(e) = /2(e) for all 
e, and the solution of Eq. ( J6 . 1 3| ) is not degenerate. 



APPENDIX D: 



He re we prove that the solution of Eqs. ( 6.13| ), ( |6.9|) , 
and (6.14) is not degenerate. Suppose the opposite is 
true, i.e. there are two functions, fi(e) and /2(e), which 
satisfy these equations. The n the function fd(e) = 
fi(e) — /2(e) satisfies Eq. ( 3.13] ) with the conditions 



def d (e) = J deef d (e) = 0, (Dl) 
f d (+oo) = / d (-oo) = 0, -1 < f d (e) < 1. (D2) 



According to Eqs. ( |6.13 ), and (|6j|), the integral 
J de\efd(£)\ is conve rgen t. This allows us to symmetrize 
the kernel a in Eq. ( |6T~3|) : a(e) = a s (e) + (1 + X)e + Ap, 
where cs(e) = [A0(e) + 6(e — Ap)] + [e — > —e] , and thus 
(Ts(e) — <r s (—e). Using the condition Eq. (Dl) we arrive 
at the new integral equation for f d , 

J de'a s (e' -e)f d (e') 

= [(l + \)(e 2 + 2T) + (Ap) 2 }f d (e). (D3) 

Next we apply Fourier transformation to both sides of 
this equation and introduce the function 



(D4) 



Here we have to be careful because, strictly speaking the 
Fourier transform of 175 (e) does not exist (this function is 
divergent at ±00). On the other hand, since the integral 
on the lhs of Eq. (|D^) is convergent, we can regularize 
the kernel as (75(e) — > <7s(e)e _r, ' e ' and later take the limit 
rj — > +0. Then for the Fourier transform of Eq. (D3) we 
find 

(1 + \)<p"(x) = [u(x) + (Ap) 2 + 2(1 + \)T]ip(x), (D5) 
u(x) = ( dee- tex a s {e) = 2[\ + cos{Apx)]/ x 2 , (D6) 



where u(x) is real, because as is an even function of 
e. Thus we have obtained a second order differential 
(Schrodinger) eq uati on for the function <p(x). We con- 
clude from Eq. (|d|) that <p(0) = (p'(0) = 0, and the 
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